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* 
Neutron Dose and Fluence Distributions in an Infinite Air Medium 


*k 
George L. Simmons and Charles Eisenhauer 


The moments method is applied to the problem of calculat- 
ing the neutron dose and fluence distributions in an infinite 
medium of air. These calculations are compared with Monte 
Carlo and Discrete Ordinates (s_) results. Simple parametric 
representations for the distributions are given which facili- 
tate the calculation of dose and flux distributions in air 


with a different density. 


Key words: Benchmark problems; dose distributions; moments 


method; neutron penetration; shielding; weapons radiation. 


I. INTRODUCTION 

The moments method for solving neutron and gamma-ray transport 
problems is ideally suited for obtaining simple, but accurate representa- 
tions of neutron dose and fluence distributions. In this note, we present 
data on neutron transport in air for three point isotropic sources of 
interest in shielding applications. These sources are: 1) fission source, 
2) 14-MeV source, and 3) a typical thermonuclear source. These results 
are compared with similar data obtained using the Monte Carlo and Dis- 
crete Ordinates (S_) techniques. We also give simple parametric repre- 
sentations for the dose and fluence which may be used to calculate these 
quantities at various distances in air of arbitrary density. 


The cross section data used in these calculations are taken from 
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the ENDF/B-I compilation for oxygen [1}' and nitrogen [2] and are the same 
as those used by Straker in Shielding Benchmark Problem 3.0 [3]. Air 
was assumed to be 79% nitrogen and 21% oxygen by volume, yielding nuclear 


densities of 4.25 x Opie angus 13 ox ar nuclei/barm-cm, respectively, 


for an air density of 1.29 grams/liter. Our calculations require that 
some prescription be given for interpolating on the basic cross section 
data. We used the specification given on the ENDF/B data file, namely, 
linear log-log interpolation on the cross sections and linear interpola- 
tion on the Legendre expansion coefficients of the elastic scattering 
angular distribution. 

Calculation of spatial distributions by the moments method involves 
two steps. First, spatial angular moments of the fluence are obtained, 
These moments are then used to reconstruct distributions at various 
depths. We used the computer code MOMENT-I [4] to calculate spatial- 
angular moments at 510 energies between 15 MeV and Tone eV. The solution 


grid used in these calculations is given in Table I. All of the source 


Table I. Solution Energy Grid Specification 


Upper Intervals 
Energy (MeV) Lethargy Number 
115/40 0.0314 35) 
350) 0.0196 200 
0.1 0.0316 60 
0.015 % 0.0769 12> 
1,0x10_44 0.1023 90 

Ox 10 


spectra used in these calculations were normalized to one neutron emitted. | 


The fission source energy distribution used has the form 


eae ; Pee , 
Figures in brackets indicate the literature references at the end of 
this paper. 
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S(E) = 0.453 Sinh(/2.29E) exp (-E/0.965). (1) 
The 14 MeV source was assumed to be represented by a histogram between 
15 and 12.2 MeV. The thermonuclear weapon spectrum is given by Straker 
[5] and was also represented in histogram form. 

In obtaining the fluence from the moments, we could perform recon- 
structions at all 510 energies. However, this would produce more data 
than can be easily managed. Also, since many of the calculational tools 
used in shielding analysis produce multigroup results, it is desirable 
to give our results in a comparable form. For these reasons, we have 


found it convenient to define multigroup moments as 


fi Ap ale ay Aa) 
oe 5 (E,-E ) | on 9 (EN dE ? (2) 
Bey in SOS Ds ees 


where E. and E. are the group boundaries for the 2 eaesey group 


j-1 
(see Table II). Consistent with this procedure, we can weight these 


multigroup moments by a dose response functions, oe and sum over all 


groups: 
2 
a) j 
D ry AE, do. oF @) 


thereby obtaining dose moments (See Table II for a set of dose response 
functions for several detectors of interest). The dose moments can then 
be used to reconstruct dose distributions. 

We have applied two different procedures in reconstructing distri- 


butions from moments: 1) function fitting and 2) plural series approxi- 
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mations [6]. The first procedure has the advantage that the representa- 
tions obtained for the distributions are relatively simple and are 
therefore easily used in point-kernel type calculations. The plural 
series technique, on the other hand, has the advantage that it enables 
us to estimate the truncation error associated with using only a finite 
number of moments in reconstructing the distribution, The function 


fitting representations make use of the form 


Shee te -E R/B, 


2 r fe) 7 Kn, 10 
4nR D(R) = Die + a, © R/B,) e (4) 
i=l 
while the plural series representation uses the product of exponentials 


and U-polynomials [6] and has the form 


-u R 


fo) Bas 


4nR7D(R) = Die (5) 


3 10 

kK k+1 k a 

+ OR) : A,/k: » a, Us R/X,)e 
ie i=1 


Here BG is the mean-free-path for neutrons at the highest energy con- 
sidered, i.e. 15 MeV, and k is a parameter which is specified for the fit 
and is usually set equal to unity. In both cases, Do is the value of 
Arr? times the unscattered flux or dose at R = Q, Do is easily calcu- 
lated from Equations (2) and (3), provided that the source energy 


distribution, S(E) is substituted for the moments, 9 49): 
II. RESULTS 


We first compare our results with existing calculations. Theoreti- 


cally, in the reconstruction of distributions from moments, an infinite 


as 


number of moments are required to obtain the correct distribution. Since 
only a finite number of moments (six to ten) are available to charac- 
terize the distribution, we estimate the truncation error due to neglect 
of remaining moments. 

The apparent convergence of the distribution, using successively 
more moments in the fit, has been used in the past to assess the trunca- 
tion error. That is to say, distributions are constructed using N 
moments and N+l moments, and the difference between the two distributions 
is examined. Generally, as more moments are used, these differences 
become smaller. However, there is no theoretical justification for 
assuming that the distributions converge to the correct value. 

Spencer [6] has developed independent procedures for estimating the 
truncation error when plural series approximations are used to represent 
the distribution, These procedures are based on a more rigorous theory 
which gives criteria guaranteeing convergence, and which gives rigorous 
bounds to the truncation error which involve the norm of the unknown 
function, We had no difficulty satisfying the convergence criteria; but 
as in Ref. 6, we rely on estimations rather than rigorous values for the 
norm and therefore also for the truncation error bounds. The truncation 
error bounds represent maximum errors, rather than probable errors. They 
are therefore a conservative estimate of the truncation error. 

For the data presented in this report, we have relied heavily upon 
the convergence of the distribution as the method for selecting the best 
fit to the moments. We have then applied Spencer's approach to determine 
rigorous error bounds on the selected fit. In Table III, we give a 


typical distribution which illustrate both the convergence and error 


6. 
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bound data for the distribution. Under the heading CONVERG, the percen- 

tage differences between an 8 and 10 moment fit is given. Note that the 

truncation error bound (BOUND) calculated as indicated in References [6] 

and [7], is larger than the convergence indicator (CONVERG). This gives 

us confidence that when convergence is good, the convergence indicator is 
a realistic estimate of the error that can be attributed to distribution 

construction, However, if this indicator is large, then we must rely 


solely upon error bound calculations to indicate the truncation error. 
A. Benchmark Problem (Air density 1.29 g/liter) 


We have calculated the multigroup fluence and Henderson dose at 
several distances for the fission and 14 MeV sources, and compared them 
with Monte Carlo calculations by Straker [3|}. In Table IV, we present 
a comparison of the multigroup fluences for the two sources for various 
distances from the source point. Generally the two methods compare very 
favorably, except in regions where the statistical error of the Monte 
Carlo results is large. Table V shows a similar comparison for the 


Henderson dose. 


B. Comparisons with Ses Calculations (Air density 1.11 g/liter) 


A large tabulation of fluences and dose has been made by Straker [5] 
for the three source types of interest, as well as for other source 
distributions. These calculations were made using the ANISN code (S,6Ps) 
for a 5000 meter sphere. Figures 1, 2, and 3 give comparisons of the 


multigroup spectra for the three sources of interest, as calculated by 


the so and moments methods. The major discrepancies are in the groups 


8. 
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Figure 1. Comparison of moments and S, (Ref. 5) calculations 
from a fission source at two depths in air. 
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12. 


Ane” x Henderson Dose (adeen-/seuree Neutron) 


R Fission Source 14 MeV Source 
(meters) Moments Monte Carlo Moments Monte Carlo 

71 4.15-09 3.98-09 5.74-09 5.97-09 
142 4.59-09 4.60-09 5.60-09 6 .03-09 
284 3.83-09 3,70-09 5.33-09 4.75-09 
380 2.93-09 2.65-09 3,72-09 3.86-09 
472 2.14-09 1.98-09 2.92-09 2.99-09 
567 1.47-09 1,32-09 2.19-09 2.20-09 
662 9.85-10 8.40-10 1.61-09 1.60-09 
755 6.49-10 5.60-10 1,16-09 1.14-09 
850 4,17-10 3.30-10 8.18-10 7.80-10 
1000 2.03-10 1,75-10 4.58-10 4.40-10 
r35 1.04-10 7.51-11 2.66-10 2.42-10 
1275 5.17-11 4.25-11 1.48-10 1.38-10 
1420 2. 84-11 2.18-11 7,.98-10 7 Spal 


Table V. Comparison between Monte Carlo and moments 
methods calculations of Henderson dose for 
a point isotropic fission and 14 MeV source 
in air (density 1.29 g/liter). 


i; 


6 and 9. At these energies there exist minima in the air cross section 
which may be difficult to represent in the multigroup procedure [5]. 
Consequently, the streaming through these minima will not be calculated 
correctly. This effect will become more pronounced as the source- 
detector separation is increased. In Table VI we show a comparison of 
the Henderson and Snyder-Neufeld dose for the three sources of interest. 
Generally the two calculations compare favorably, with the moments 
results tending to be higher than the so results. We feel that this is 
due to the increased streaming that is predicted by the moment calcula- 
tions, as well as to the leakage effect that is associated with the 


finite sphere used in the so calculations. 


C. Simple Representations for Fluence and Dose Distribution. 


For point kernel calculations, it is useful to have a simple set of 
parameters which can be used to represent the dose and fluence for the 
three sources. Since we are interested in simple representations of 
these distributions, the results of the function fitting calculations 
are used, We note in passing that function fitting and plural series 
techniques give the same distribution over the range of interest. In 
making these tabulations, we restrict ourselves to representations of 
the distributions by 11 or fewer parameters. In Tables VII, VIII, and 
IX, we give coefficients for the dose quantities in Table II. The form 


used in calculating the distributions from the coefficients is 


5 
ER 5 R/B. 
4mR°D(R) = Die ° +) a, 0 R/B,)e ° * (rads-em*) (6) 


i=l 


14, 


AnR2 X Henderson Dose (Rads-cm“/Source Neutron) 
Thermonuc lear 


R Fission Source 14 MeV Source Weapon Source 
(Meters) Moments s, Moments so Moments so 
150 4.50-09 4.49-09 5.69-09 6.07-09 2.41-09 2.49-09 
300 4.05-09 3.93-09 4.88-09 5 .06-09 1.95-09 2.03-09 
600 1.81-09 1.76-09 2.59-09 2.60-09 8.55-10 8.85-10 
900 5.92-10 5.82-10 1.09-09 1.08-09 3.08-10 3.13-10 
1200 1.72-10 1.68-10 4.04-10 3.94-10 1.02-10 1.01-10 
1800 129-11 en G1: 4.48-11 4.24-11 OS U2 eo. 2 = 1 2 
2400 DESO US NM RES} 4,32-12 3.85-12 
3000 6.80-14 4,.27-14 JO ES Sel B/E) 
3600 5.08-15 2.52-15 3.34-14 2-46-14 
4200 3.84-16 1.49-16 DMS 1.84-15 
4800 2.92-17 8.43-18 2.32-16 1.26-16 
hap? X Snyder-Neufeld Dose (Rads-cm“/Source Neutron) 
Thermonuc lear 
R Fission Source 14 MeV Source Weapon Source 
(Meters) Moments so Moments so Moments so 
150 6.80-09 7.40-09 8.40-09 8.25-09 5.47-09 5.49-09 
300 7.08-09 TALS SON) 7.01-09 7.22-09 4.42-09 4.55-09 
600 3.82-09 3.62-09 4.05-09 4.05-09 1.76-09 1.85-09 
900 1.30-09 1.28-09 1.84-09 1.79-09 6.15-10 6.39-10 
1200 3.86-10 3.82-10 7.10-10 6.88-10 2.01-10 2.02-10 
1800 2.89-11 2./70-11 8.24-11 7.8/-11 OS nel 1.84-11 
2400 2.02-12 1.66-12 8.11-12 7.39-12 
3000 1.41-13 9.78-14 7.35-13 6.20-13 
3600 1.03-14 5.68-15 6.36-14 4.86-14 
4200 8.10-16 3.31-16 5.35-14 3.65-14 
4800 6.84-17 1.83-17 4 .40-16 2.48-16 
Table VI. Comparison between S, calculations and moments method calcula- 


tion of Henderson and Snyder-Neufeld dose for a point isotropic 
fission, 14 MeV, and thermonuclear weapon source in air (den- 
Sivas bmeienber)). S, data taken from Reference 5. 


ie 


CFONEE CFUNFEE 


Oo fFuUNEF FE 


HENDERSON (RAD) 


Do 72671191009 
ay 
=4e515538-09 
8e190579=09 
5°651962=-09 
2e513979-09 


K= 1 


1.67617 


Do : MID=PHANTOM (RAD) 


8-/76196-10 
az 
4eU18550-09 
4e205045=U9 
-9°U05U 73-09 
12109593=10 

0-UQ0U00U0 


NON=10Ne 
Do :5¢/10731-09 
az 


9-0U3585=-10 
=1¢535754=-08 
1¢402255-08 
1¢073952-08 
20¢475567-09 


4nR X 


Table VII. 


KS 1 
Bi 
°49140 
1.47419 
66646 
1285908 
1.00000 


SILICON KERMA 


K= 1 


1644839 
1.69788 


FISSION SOURCE IN AIR 


SNYDER=NEUFELD (RAD) 


3-843843-09 K= 1 

Og Bi 
3-032895-08 .48881 
22651148-08 1.46643 
-3+700481-08 63225 
60949531-10 1.85179 
0.000000 1.00000 


CONCRETE KERMA (ERG/G) 
30453603-08 K= 
aq Bi 
9-839580-09 
-5-994808-08 
32863627-08 1.39454 
12162011-07 1235447 
3+363860-08 167387 


IONe SILICON KERMA 
42.329498-09 K= 1 
ad By 
723547080=-10 
-4.782786-09 
1.069465=08 
22156281=-09 
3-044102=09 


°88542 
1237297 
1.00379 
1.66981 


5 


i= 


a point isotropic fission source. 
have units of (dose units-cm2 per source neutron). 
133.059 meters for an air density of 1.11 g/liter. 
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TISSUE KERMA 
2°752741=07 
aj 
9-387321=08 
=5-611913=07 
8+124505=-06 
-6-808031=06 
3-2164664=07 


(ERG/G) 
K= 


»90590 
1640472 
1641258 
1.66805 


AIR KERMA (ERG/G) 


4.324899=08 
Oy 
9-665720-08 
1-955809=07 
=-1+497707=07 
6-125755=-09 
0.000000 


°49046 
1247138 
©77937 
1.84231 
1.U0000 


SNYDER@=NEUFELD (REM) 


3210515-08 
ag 
-32173191-07 
4.050881=-07 
5+«584109=08 
0.000000 
0.000000 


k 
Dose(R) = D,exp(-E4R) +). a, @R/B,) exp(-2R/B, ) 
a 


K= & 
Bi 


12135018 
1.25199 
1.04183 
1.00000 
1.00000 


Coefficients for calculating dose distributions in air for 
The dose distributions 


XS is 


uUFrOonRe uFfonrr 


uFfOoNrY 


Doz 


NON*LONe 
D,: 


Table VIII. 


HENDERSON 


9°40UU019=-09 
ai 
-60645252-10 
92663185709 
6¢827104-09 
-1-¢4276/7-08 
9e676155-09 


(RAD) 
Ka ad 
Bi 
225261 
1.05069 
1.62923 
1214732 
1.72521 


M1D=-PHANTOM (RAD) 


4e9UU0017-09 
CsI 

4. /43589=-10 

le /21504=-09 

52092554=09 

-2e612751-09 

4e960820-09 


7e5U0026-09 
al 
2¢/54358-08 
-1+5686010-08 
12971405=-09 
-5-8604663-12 
0-V00UU0 


1.71004 


SILICON KERMA 


K= 1 
Bi 
1.63360 
1.40440 
1.86455 
2655862 
1.00000 


14 MeV SOURCE IN AIR 


SNYDER=NEUFELD 
7.000024-09 
ay By 
-8-955006-09 
7.164796-09 
4898157-08 1634196 
2.860304-08 170389 
-6+504560-08 1631416 


(RAD) 


CONCRETE KERMA (ERG/G) 
1.580005=-07 
ay Bi 
-1.517908-08 
6+707914=08 
8-164871-08 1.54762 
-9.334668-08 1.16231 
1.620361=-07 1.71573 


IONe SILICON KERMA 
8-600030=-08 
ay By 
—6e418379-09 
7 »669295-09 
4e430919-08 1.35636 
32612633-08 1259763 
72979537-09 1474196 
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TISSUE KERMA 
6¢360022=-07 
Ss 
-7+657209-08 
1-172000-00 
8e068461-07 
-1-694500-06 
1-003330-06 


AIR KERMA 


2-660009-07 
eet 
-2-829879-08 
72271133-08 
-1.071440-07 
22252018-07 


(ERG/G) 
K= 


1.72744 


(ERG/G) 
K= 1 


1.50521 
eS ier 
1.71239 


SNYDER@=NEUFELD (REM) 


5-640019-08 
es 
-9-282823-09 
4.560220-07 
2e022771-07 
=5-393797=-07 
1.088691-08 


= k 
Anke X Dose (R) = Dexp(-Z,R) +) aj ZoR/B,) exp (-2,,R/B, ) 
i=1 


1.61020 


Coefficients for calculating dose distributions in air for 


a point isotropic 14 MeV source. The dose distributions 


have units of (dose units-cm2 per source neutron). eer: 
133.059 meters for an air density of 1.11 g/liter. 


izle 


OCOfonrre yftonrr 


vu FUNPe 


THERMONUCLEAR SOURCE IN AIR 


HENDERSON (RAD) 
Do: 1¢598693-U9 K= 
ot Bi 

262355U79-10 53182 
22229882=09 1455445 
20699256-09 1.69711 
2¢951245-10 1275437 
UeOQU0UU0 1.00000 


MID=PHANTOM (RAD) 
Do: &¢9U5867-10 K= 1 
aj Bi 
-72557294-10 21165 
12531144=-09 .88030 
121600956-09 1.36502 
=-8-230U29-10 .99819 
16135242-09 1.79155 


NON=IONe SILICON KERMA 
Doi 2e104689-09 K= 1 
ag Bi 
—-2-849184=-09 .4A294 
1-5/8367°U9 1.44882 
2e592974=09) 43650 
167718536709 1-28652 
53285554609 1.70903 


SNYDER=NEUFELD (RAD) 
20462533-09 K= 1 
ag Bi 
5-087558=-09 1237554 
-8.094390-09 .47979 
96294106-09 .63088 
9-040126-11 1.60195 
5-399257-09 1.70926 


CONCRETE KERMA (ERG/G) 
26863221-08 K= 1 
az Bi 
4.211789-09 .61293 
36432219-08 1235017 
3-735728-08 1.69671 
2.190751-09 1-77446 
0.000000 1.00000 


IONe SILICON KERMA 
12139916-08 K=_1 
az Bi 
-1.102444=-10 1.35895 
=-32306599-12 .40243 
32814736-09 1.17663 
7e002446=-09 1445429 
4.738565-09 1.68985 


5 k 


TISSUE KERMA 
1-772189-07 
SL 
-3-418369=-028 
5-063562=-U08 
2°976032=-07 
3-183777-07 
0.000000 


AIR KERMA (ERG/G) 


4.131295-08 
aj 


42733615-08 
-7-026025-08 
7-407690-08 
4.2939778-08 
0.000000 


SNYDER=NEUFELD (KEM) 


1.842464-08 
aj 


-1.560135-08 
2¢265075=-08 
1-859109=-08 
1.240604=-08 
3e¢849104-08 


AnRe X Dose(R) = Doexp (-29R) 4) a; @oR/B; ) exp (-2R/B; ) 
i=1 


Table IX. Coefficients for calculating dose distributions in air for a 
PTI thermonuclear weapon source. 
units of (dose units-cm2 per source neutron). % is 133.059 


The dose distributions have 


meters for an air density of 1.11 g/liter. 
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(ExG/G) | 
K= 1 
Bi 
1.41642 | 
-ou201 
1.373511 
1.70691 | 
1.00000 | 


Kent 
Bi 
1.39377 
57382 
° 82382 
1.70684 
1.00000 


1.71016 


Here 1/d , = 133.905 meters which corresponds to an energy of 15 MeV and 
an air density of 1.11 grams/liter. These same coefficients may be 
applied to calculations at different air densities by using simple scaling 
procedures based on equivalent mean free paths, In Tables X, XI, and 
XII, the distributions obtained by using the coefficients in Tables VII, 
VIII, and IX are given. 

In Tables XIII, XIV, and XV, the coefficients which represent the 
multigroup spectra for each source are given, The form for these dis- 


s : Zien 2 
tributions for the j energy group is 


7 
LR ie > R/8. My 


4nR7g7 (R) = 9 oy +’ a, @ R/B,)e . * (neutrons/MeV/source neutron) 
i=1 

Tables XVI, XVII, and XVIII give the distributions obtained using these 
coefficients. Note that for the thermonuclear weapon source, at very 

low energies and close to the source, the distributions are negative. 

This indicates that the coefficients for these particular energies should 
be used with caution at distances less than one mean-free-path. Indeed 
in this particular region, it is extremely difficult to fit the moments 


and obtain suitable distributions over the entire range. 


D. Angular Distributions. 


Angular distributions for dose or fluence can also be calculated by 
the moments method, Generally the procedure that is used for point 
Sources is to calculate spatial distributions for Legendre harmonic 
coefficients py (E, 2) for values of £ up to about 8. In theory,-the 


angular distributions at any distance can then be calculated by perform- 
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ing a Legendre sum of the harmonic coefficients for the appropriate 
values of E and z. In practice, the number of coefficients available is 
insufficient - particularly for high energy neutrons - and spurious 
oscillations or even negative values occur. These oscillations can be 
reduced considerably by using an extrapolation technique devised by 
Spencer [6] for estimating harmonic coefficients for higher ¢ values. 
Although this technique does not eliminate oscillations, the remaining 
peculiarities are very often limited to negative values of cosQ, where 
the angular distribution has relatively low values and is therefore of 
little practical interest. The angular distribution near the source, 
particularly for point sources, tend to exhibit trends associated with 
the single scattered angular distribution. Thus near the source, one 
expects that the angular distribution would exhibit oscillations similar 
to those of the scattering kernel. 

In Figure 4 we show a comparison between moments and So calcula- 
tions of the Snyder-Neufeld dose angular distribution for a fission 
source. These comparisons indicate that the moment results tend to be 
more peaked forward than the Ss results. Consequently more neutrons 
would transport further from the source. This phenomena also contributes 
to the higher dose values at large distances that is predicted by the 


moments calculations. 
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SNYDER NEUFELD RAD DOSE ANGULAR DISTRIBUTION 
POINT ISOTROPIC FISSION SOURCE IN AIR 


—— MOMENTS 
= EST 


Air R2 D(R,COS @) [RAD x cm*/ (STERADIAN-SOURCE NEUTRON)] 


COS @ 


Figure 4. Comparison of moments and S, (Ref. 5) calculations of the 
angular distribution of Snyder-Neufeld dose in rads at four 
depths in air. 
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